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Abstract 

This text presents the proofs of the technical facts underlying theoretical justification of the convergence and performance of the novel 
algorithm for reactive navigation of differential drive wheeled robots in dynamic uncertain environments. The algorithm restricts neither the 
natures nor the motions of the obstacles, they need not be rigid but conversely may deform. It does not consume data about the velocities, 
shapes, sizes, or orientations of the obstacles, and does not need a map of the environment or recognition of individual obstacles. The 
only information about the scene is the current distance to the nearest obstacle. 



> 

in 

t> 

^t 

o 
m 



% 



1 Introduction 

A key component of safe navigation is avoidance of colli- 
sions with en-route obstacles. This problem has attracted an 
enormous attention in robotics research. The relevant algo- 
rithms can be generally classified into global and local path 
planners [1]. 

Global planners typically build a more or less comprehen- 
sive model of the environment to find the best complete 
trajectory [2]. By and large, they are computationally ex- 
pensive; NP-hardness was established for even the simplest 
problems of dynamic motion planning [3]. This seriously 
troubles their real-time implementation. Data incomplete- 
ness and erroneousness, typical for onboard perception, may 
cause a noticeable deterioration in the overall performance 
of global planners. 

Conversely, local planners iteratively re-compute a short- 
horizon path based on sensory data about a close vicinity of 
the robot. A short computation time typical for these plan- 
ners creates a potential for their use in real-time guidance 
systems. However most of the related techniques either in 
fact treat the obstacles as static or assume a deterministic 
knowledge about the obstacle velocity and a moderate rate 
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of its change; the first group is exemplified by the dynamic 
window [4,5], the curvature velocity [6], and the lane cur- 
vature [7] approaches, whereas velocity obstacles [8], colli- 
sion cones [9], or inevitable collision states [10,1 1] methods 
provide examples of the second kind. 

In the marginal case where the planning horizon concen- 
trates into a point, local planner acts as a reactive con- 
troller: it provides a reflex-like control response to the cur- 
rent observation. Examples include artificial potential ap- 
proach [12,13] and kinematic control based on polar coor- 
dinates and Lyapunov-like analysis [14]. In this area, the 
obstacles were mostly viewed as rigid bodies of the sim- 
plest shapes [13,14,12,15], the sensory data were assumed 
to be enough to provide access to the locations of the ob- 
stacle characteristic points concerned with its global geom- 
etry (e.g., the disc center [13,14] or angularly most distant 
polygon vertex [15]) and to the full velocity of the obstacle 
[13,15,14]. Furthermore, rigorous justification of the global 
convergence of the proposed algorithm was rarely encoun- 
tered. 

In this paper, the problem of reactive navigation is addressed 
for a standard differential drive wheeled robot. The objective 
is to drive the robot to the assigned point-wise target through 
the obstacle-free part of the scene. The robot may have a 
limited sensor range, may not distinguish between various 
points of the obstacle, and so may be unable to estimate many 
of its parameters, like size, center, edge, full velocity, etc. It 
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has access only to the distance to the nearest obstacle point 
within the sensor range, the time derivative of this reading, 
and the angle-of-sight at the target. Unlike the overwhelming 
majority of the previous works, the obstacles are not assumed 
to be rigid or even solid: they are continuums of arbitrary 
and time-varying shapes undergoing general tranformatons, 
including rotations and deformations. 

The extended introduction and discussion of the proposed 
control algorithm are given in a paper submitted by the au- 
thors to the IEEE Transactions on Control Systems Technol- 
ogy. This algorithm was extensively tested via both computer 
simulations and real-world experiments with both research 
wheeled robots and an intelligent wheelchair. The results of 
both simulation tests and experiments are also presented in 
that paper. This text basically contains the proofs of the tech- 
nical facts underlying theoretical justification of the conver- 
gence and performance of the proposed algorithm. These 
proofs were not included in that paper due to the length lim- 
itations. To make the current text logically consistent and 
self-contained, we supply the reader with relevant notations 
and complete statements of both the problem and theorems 
to be proved in this text. 

2 System Model and Problem Statement 

We consider a planar differential drive wheeled robot 
(DDWR) with two independently actuated driving wheels 
mounted on the same axle and castor wheels. DDWR is 
controlled by the angular velocities uji and cu r of the left 
and right driving wheels, respectively. These velocities are 
limited by a common constant 51 > \u>i\, \u> r \; the driving 
wheels roll without sliding. 




Fig. 1. Differential drive wheeled robot 



The kinematics of such vehicles are classically described by 
the following equations: 



R w is the radius of the driving wheels, 2L is the length of 
the axle, and uji = u>i(t) G [— fi, fi], i = I, r. To simplify the 
matters, we treat v and u as control variables. They uniquely 
determine the rotational velocities uj r = (v + Lu)/R w ,uji = 
(v — Lu)/R w and obey the bound: 



L\u\ < V := R w n. 



(2) 



It follows that for given v G (— V, V), the turning radius of 
DDWR is no less than 



R 



L\v\ 



V 



(3) 



The workspace of DDWR contains a moving obstacle 
D(t) G K 2 , which need not be a rigid body but con- 
versely, may undergo arbitrary motions, including defor- 
mations. Its position and shape are not known in advance; 
the only available information is the current distance 
d(t) := disto(t)[)~(£)] t0 me obstacle and the rate d(t) at 
which this reading evolves over time. Here r := [x, y] T is 
the vector of the coordinates of DDWR and 



dist£)(r) :~ min \\r — r 

r'GD 



(4) 



where || • |j denotes the standard Euclidean vector norm and 
the minimum is achieved if D is closed. 

Finally, there is a steady point-wise target T in the plane 
and DDWR has constantl access to the heading h(t) towards 
the target. The objective is to guide DDWR through the 
obstacle-free part of the plane and reach the target T at a 
certain time tf > 0: 



'(*/) 



r(t)#D(t) Vte[0,t f ]. 



Moreover, the distance to the obstacle should constantly ex- 
ceed a given safety margin d sa f e > 0: 



dfatc W [r(t)} > d S afe Wte[0,t f ] 



(5) 



x = v cos ( 



2 ' 



;(0) 



XQ 



y = vsin9, u=^ r ^ L , y(0) = y 
= u, Vi = R w uJi, 9(0) = 9o 



(1) 



Here x, y and 8 are introduced in Fig. 1 and stand for the po- 
sition and orientation of DDWR, v is its longitudinal speed, 



3 Autonomous Navigation Algorithm 



Now we present a summary of the proposed reactive nav- 
igation algorithm. The algorithm combines obstacle avoid- 
ance behavior, which is activated in a close proximity of the 
obstacle, with motions towards the target in a straight line 
when there is no threat of collision. 



We employ the following obstacle avoidance strategy: 



u(t) = ^^-sgn{d(t) + x[d(t)-d }} 



V-v(t) 
L 

v(t) = T[d(t)] 
X(z) : 



, where (6) 



72: if \z\ < 5 

iJ*sgn(z) if \z\ > S 



(«♦ := 7<5) (7) 



is the linear function with saturation, the smooth function 
T(-) : [0,oo) -)• (0,V) determines the speed u of DDWR on 
the basis of the current distance D to the obstacle, whereas 
7 > 0, 8 > 0, and do > d sa f e are the controller parameters. 
Here do has the sense of the desired distance to the obsta- 
cle when bypassing it. The function T(-) varies between 
two speeds vq and v CI , i.e., T(d) = v Q Vd < dJ,T(d) = 
«cr Vd > d cr > do . The speed v e (0, V) is used when by- 
passing obstacles, so do < dj ; the larger v a - > vq 'cruise' 
speed is employed where there is no collision threat. We 
assume that v a < V to leave the vehicle a certain level of 
maneuverability in the "cruise" regime. 

The proposed obstacle avoidance strategy belongs to the 
class of sliding mode control algorithms; see e.g. [16]. The 
intuition behind this strategy is that in the sliding mode, the 
equation d+\(d — do) = of the sliding surface is satisfied, 
according to which the vehicle is steered towards the desired 
distance do to the obstacle: d — > do. For this to take effect, 
the sliding mode maneuver should be at least realistic. Since 
the derivative d does not exceed \d\ < v r the relative speed 
v r of the vehicle with respect to the obstacle, this means that 
in (7), the saturation level v* should not exceed this speed. 
This can be achieved by proper tuning of the controller pa- 
rameters 7 and S based on available estimates of the speed 
of the obstacle. If initially the vehicle is not on the sliding 
surface, the control law (6) quickly drives it to this surface 
after a short initial turn, see Section 6 for details. So sliding 
motion is the main part of the obstacle avoidance maneuver. 





Fig. 2. (a) Sliding motion towards the equidistant curve; (b) Motion 
along the equidistant curve 

In more details, this motion looks as follows. The equation 
d + x(d — do) = means that the angle a between the 



relative velocity of the vehicle v r and the line of sight at 
the nearest point of the obstacle equals a — arccos x ^ - — . 
Hence the angle a is obtuse for d < do and acute for d > do, 
and so the vehicle is driven towards the desired distance 
do to the obstacle; see Fig. 2(a). In doing so, a is kept 
constant a = arccos[— sgn%(d— do)] in the saturation zone 
d — do I > S. As d leaves this zone and approaches do, the 
angle goes to J. In the limit where d = do, the vehicle is 
oriented parallel to the obstacle boundary a = ir/2, which 
means traveling along the do -equidistant curve; see Fig. 2(b). 

The control law (6) is activated in a close proximity of the 
obstacle. Otherwise the vehicle is driven towards the target 
in a straight line: 



u(t) = 0, v(t) = T[d(i)]. 



(8) 



Switching (8) M- (6) occurs when the distance to the nearest 
obstacle does not exceed a given threshold d av G (dojdj] 
and d + x(d — do) < 0; switching (6) n- (8) occurs when 
the vehicle is headed for the target and d + x(d — do) > 0. 

The second relation from (8) does not mean that the vehicle 
should have constant access to the distance d(t): it suffices 
that this distance can be measured only in a vicinity of the 
obstacle d(t) < C, where C > d cr . 

Practically the discrepancy d cr — d av and the profile of T ( ■ ) on 
[d av , d cr ] are chosen with regard to the acceleration capacity 
of the wheelchair. 



4 Main assumptions 

To describe the properties of the moving continuum D(t)P1 
we introduce a reference configuration D* C M 2 of the body 
and the configuration map <&(•,£) : K 2 — > M 2 that transforms 
D* into the current configuration D(t) = $[£>*, t). We limit 
ourselves by only few conventions typical for the continuum 
mechanics and listed in the following. 

Assumption 1 The reference domain D* is compact and 
has a smooth boundary dD*. The configuration map $(•, t) 
is defined on an open neighborhood of -D* and is smooth, 
the determinant of its Jacobian matrix is everywhere positive 
dct & r > 0, and $(ri, t) ^ $(t"2, i) whenever r\ 7^ r^. 

Hence the boundary dD(t) is smooth at any time. 

To proceed, we need some notations, which partly refer to 
the Eulerian formalism in description of the body motion: 

• < I >_1 (-, i) — the inverse map; 

• V(r,t) := $t[$ _1 (r,t)] — the velocity vector-field; 



* Discussion of the properties of the moving continuum D{i) 
basically follows that in [17]. 



A(r,t):~- 



$" t [$ 1 (r, t)] — the acceleration field; 
the spatial velocity gradient tensor; 



v;[r 7 t] + v;[r : t] T 



the strain rate tensor; 



uj(r,t) — the vorticity, or angular velocity of the 



o u> 

-u) 



V' -V 



rigid-body-rotation i.e., 
(•; •) — the standard inner product; 
f(r,t),N(r,t) — the Frenet frame of dD(t) at r G 
dD(t); the domain D(i) is to the left to the tangent vector 
T, the normal vector N is directed inwards D(t); 
Wr,t(r,t) := (W(r,t);f(r,t)) — the tangential com- 
ponent of the vector-field W at r G dD(t) at time t; 
WN,t(r, t) — the normal component; 
re(r, t) — the signed curvature of dD(t) at the point r; 

a[r, t] := (£T ; N) + u = (VJ.T; N) — the normal com- 
ponent of the velocity rate-of-change under an infinitesi- 
mally small shift along the boundary dD(t); 
r(t) — the current position of the vehicle; 
7%(£) — the point of dD(t) nearest to r(i). 



The signed curvature k assumes positive and negative values 
on the convex and concave parts of dD(t), respectively. 

The next assumption is typically fulfilled in real world, 
where physical quantities take bounded values. 

Assumption 2 The scalars V/v[r,£], Vr[r,t], A N [r,t], 
a[r,t], K,[r,t] remain bounded as r G dD(t) and t — > oo. 

As t runs within any finite time horizon, these quantities do 
remain bounded over r G dD(t) since all of them, along 
with $(r,i), continuously depend on r and t, and the ref- 
erence domain D* is compact by Assumption 1 . 

In order that the vehicle be capable of bypassing the obsta- 
cle at the desired distance do, the do -equidistant curve to be 
traced should not be sharply contorted and move too fast. 
The first and second parts of these requirements are fleshed 
out by the following Assumption 3 and Lemma 1. The for- 
mer is shown to be nearly unavoidable even in the case of 
the static D(t) = D domain [19], whereas the proof of the 
latter can be found in [17]. 

Assumption 3 For any point r on the do- equidistant curve 
of D(t), the distance from r to D(t) is furnished by a single 
point r» G D(t), and 1 + do^fr*, £) > 0. Moreover, 



lim inf (1 + d K[r*,i\) > 0. 

£->oor„e&D(t) 



(9) 



Lemma 1 Let the vehicle travel so that it constantly over- 
takes the nearest boundary point of the obstacle and d(t) = 
do, v(t) = v. For any t, the parameters of the obstacle mo- 



tion at the point [7% (£),£] satisfy the following relations: 

\A\L 



V^Vt 



+ v<V, 



\V N \ < v, 



A = A(r*,t, d , v) := A N + 



u 2 -V%>±(V T + d a) , 



2cr£ + k£ 2 - d Q o- 2 
1 + ndo 



and £:=-V T ±\/v 2 -V$. (10) 



Here the sign + is taken if the vehicle moves so that the 
obstacle is to the left, and — is taken otherwise. 

The assumption v(t) = v follows from d(i) = do if me 
speed is a function of the distance d(t), like in (6) and (8). 
Thus conditions (10) must be satisfied, otherwise the control 
objective cannot be achieved. Since the nearest boundary 
point r*(i) is not known in advance, it is reasonable to 
require that (10) holds for all boundary points at any time. 
The next assumption enhances a bit this by substituting the 
uniformly strict inequality in place of the non-strict one. 

Assumption 4 There exist X v ,X a G (0,1), £„ > 0, v a G 
(0, V) such that for v :— vq, the following inequalities hold: 



\Vn\ < \ v v, 



\A\L 



V^^V, 



+ v<X a V, 



v 2 -Vl > \VT + doa\+e. L 



(ID 
(12) 



at any time t, point r» G dD(t), and with both signs ±. 

For any finite time horizon, such A„, A , e v exist if and only 
if the inequalities in (11) are strict. Assumption 4 ensures 
that the strict inequalities do not degrade as t — > 00. This 
assumption clearly holds for motionless obstacles (An = 
0, V T = V N = 0, a = 0, £ = ±v) with v = v « 0. 

Inequalities (11) allow not only to maintain the distance to 
the obstacle d = but also to both decrease d < and in- 
crease d > it via proper manipulation of the control u, 
thus making the output d locally controllable. As may be 
shown, the above assumptions are enough for local stabil- 
ity of DDWR in a vicinity of the desired equidistant trajec- 
tory. To establish non-local convergence to this trajectory, 
we need to ensure controllability during the transient. To 
this end, we define the launching motion as that with the 
speed vq and maximal turning rate u = —(V — vq)/L start- 
ing at a time t* from a position at the distance d av from the 
obstacle D(t*). This motion is over a circle of the radius 
voL/(V — vq) and is nothing but the motion that is typically 
observed at the initial stage of obstacle avoidance maneuver. 

Our final assumption encompasses the previous nearly un- 
avoidable assumptions. 




Fig. 3. Launching motion 

Assumption 5 There exist v € (0, V), A„, A a G (0, l),e v > 
0, d- < d+ such that the following claims hold: 

(1) d sa f e < d- < do < d + and (9) — (12) are true for 
v :— vo, do := d and any d G [d_, d+]; 

(2) During the first one and a half turn of any launching 
motion the distance d(t) from DDWR to D(t) lies in 
the interval [d_,d+]; 

(3) In this motion, the total rotation angle a of r(t) — 
r*(t) does not exceed (i* — T turn ){V — vq)/L at some 
time U G [Ttum, Ti.stum], where T tum and n.stum are the 
times when one full turn and 1.5 turns are completed, 
respectively. 



where £, := -V T ± \Jv 2 - (V N + z) 2 (14) 

is continuous in z w uniformly over r» G dD(t),t > 
0, d G [d_, d_|_] and for z := 0, equals the left-hand side of 
the second inequality from (11). Hence there exists z* > 
such that for above r*,t,d, 

n(r*,t,d,z)<(\ a +r) a )V Vz G [-z*,**]. (15) 



Now we are in a position to state the main theoretical result. 

Theorem 1 Let Assumptions 1, 2, and 5 be satisfied, the 
parameter vo of the controller be taken from Assumption 5, 
and the other parameters be chosen so that v cr G (0, V), 



v* = jS < mm{r) v v , z t }, 

(A + f] a ) + 



7-Lw* 



v {V - I'oV 1 " (K+Vv)' 



<1. (16) 



Suppose that after the distance from DDWR to D{t) reduces 
to the dangerous level d av , the avoidance control law (6) is 
activated and then drives the vehicle. Then the distance to 
the obstacle constantly exceeds the safety margin d sa f e and 
moreover, lies within the interval [d_ , d + ] from Assumption 5 

and goes to the desired value d(t) > do-, d(t) — > 0. 

Since some time instant, DDWR constantly overtakes the 
nearest point of the obstacle boundary. 



Here 3) holds with t* = r turn for steady obstacles since then 
a = 0. For moving and deforming obstacles, 3) holds with 
t* = ti.5 mm if e-g-, the obstacle D(t) and WC remain in 
disjoint steady half-planes and so a < n; see Fig. 3(a). 

Since in (10), A is the (first degree)/( first degree) rational 
function in do, simple calculus show that it suffices to check 
the second inequality from (11) only for d = d^,d + . Simi- 
larly (9) holds for all d G [d_, d+] if and only if 

lim inf (1 + d+nlr, £]) > 0. 

t^oor£dD(t):K[r,t]<0 



The proof of this theorem is given in Section 6. 

Chosen \ Vl r\ v ,v = v,a and 7, the parameter requirements 
(16) can be always satisfied by picking S small enough. 

Remark 1 Theorem 1 remains true for the navigation law 
(6) with the reversed sign 



»(*) = -^F 1 ■ sgn{d(t) + X [d(t) - do]} 



v(t) = T[d(t)] 



(17) 



5 Tuning controller parameters and the main results 

In this section, we consider the parameters A„, A a < 1, d_ < 
d+, vo taken from Assumption 5. To tune the controller 
parameters, we first pick two reals t] v > 0, i] a > so that 

A„+?fo<l, \ a + Va< 1 - (13) 

For any choice of the sign in ± and v := Vo, the function 



£l(r*,t,d,z) := L 



A N + 


2at+mC- 


dcr 2 


1+nd 




v^ 2 - 


~{V N + 


zf 



V, 



provided that in the definition of the launching motion, the 
turning rate u = (V — Vq)/L is put in place of u = — (V — 

v )/L. 

Target reaching property requires much more sophisticated 
analysis even in the case of steady obstacles. In general, the 
control law (6) does not ensure that dealing with a given 
obstacle will terminate [20]. However, this may hold only for 
complex maze-like obstacles [20]. Even in this case, there is 
a remedy in the form of a random choice among the control 
laws (6) and (17) whenever a new avoidance maneuver is 
commenced [20]. Unavoidable target reaching in the scene 
with steady convex obstacles was theoretically demonstrated 
in [19] for a unicycle driven by a similar control law. 



6 Proofs of theoretical results 

From now on, the Frenet frame [T(r,i),N(r,t)] and 
the variables attributed to the motion of the domain, like 
V(r,i),A(r,i), etc., are considered only for r := /•*(£). 
With a slight abuse of notations, the resultant argument 
[?%(£), t] is replaced by t and dropped if t is clear from the 
context. 



and two-side repelling in the sub-domain 

C + V T <0 (24) 

if (16) holds. On this surface within the above domain, 



\£ + V T \ >v^/l-(X v +r, v ) 2 >0. (25) 



Lemma 2 ([17]) The velocity v of the robot has the form 

v = [s + dp]f+V -dN, (18) 

where s is the speed of the relative motion of r* along the 
boundary dD, i,e, r* — V = sT, and 



/i := (v^f , N\ + k's = Uf, n\ +lo + k's. 

Lemma 3 ([17]) The following relations hold 

d = £/i - u [£ + V T ] + A N 
+s(v;f;N)-^(V N -d) 



(19) 



= -u [£ + Vt] + A N 

v 



2cr£ + < 2 - da 2 



1 + nd 



(V N - d); 

Z-do 
' T+~Kd' 



(20) 

(21) 
(22) 



where s is the speed of the relative motion of r* along the 
boundary dD, i,e, r* — V = sT. 

Proof of Lemma 1. Since d(t) = do =>■ d(t) = Q,d(t) = 
0, v(t) = F'[d(t)]d(t) = 0, (18) yield that v N = V N ,v T 



£ + Vt, where v T + ■ 



'N 



< V 2 , which implies the first 



and last relations from (10). Similarly we have by (21), 



u = ±- 



Vv^Vi 



A 



N 



2a£_ + < 2 - da 2 
1 + nd 



The proof is completed by noting that |u| < ^jt ^ ^' anc ' 
overtaking means that s > if the vehicle has the domain 
to the left, and that s < otherwise. 

From now on, the assumptions of Theorem 1 are assumed 
to hold and we consider the wheelchair driven by (6). 

Lemma 4 Within the domain d 6 [d—,d+], the surface 
S := d + x(d ~ do) = is sliding in the sub-domain 



C + V T >0 



(23) 



PROOF. Whenever S := d - v = 0, where v := ~x(d - 
do), we have by (7), 

Ml = M < «*; H < l\d\ = l\x(d - d )\ < 7 w*. (26) 
Due to (18), 

{£ + V T ? + {V N -df=v 2 . (27) 

and so 



|£ + V T \ = y/v* - (V N - d) 2 > Vv 2 -(\Vn\+v*) 2 



(a) 



>tVl- ^+7^)2^(25), (28) 

where (a) holds due to (11) and (16). Furthermore, 

S = —[d + x{d- do)] = 'd-i> 

(2D , t , T/ , , . , 2fx£ + < 2 - da 2 . 

= -u £ + V T \ + A N H — — v, 

1 + nd 



where £ = — Vt ± J v 2 — {Vn ~ d) 2 due to the starting 
argument from (28). Let B(-) stand for the expression within 
| • | in (14) and u := (V — vq)/L. Then we see that 



S=[i + V T ] 
Here due to (15), 

B(r*,t,d,d)- 



-u + B(r*,t,d,d) 



Z + V T 



(29) 



z + v T 

(25) _ 

< (A„ +rj a )u + 



< (A a + T]a)u + 

JV* 

vy/l - (\ v +rj v ) : 



7«* 



\S + Vt\ 

(16) 

: < U (30) 



So the signs taken by S for u = u and u = —u, respectively, 

are opposite. In the case (23), the sign is opposite to sgnit = 
sgnS*; in the case (24) the signs are equal. This implies the 
conclusion of the lemma. • 

Lemma 5 If the equation d + \{d — do) =0 becomes true 
at some time t when d £ [d_ , d + ], then monotonically and 
exponentially fast d — y do , d — > as t — > cxd. Furthermore, 
d(t) € [d-,d+] and (23) holds for all t > to- 



PROOF. Lemma 4 guarantees that first, £ + Vt > at 
t = to and second, this inequality is still valid and sliding 
motion occurs while d £ [d_, d+]. During this motion, y = 
~x(y) f° r V '■= d — do, where x{v) ' V > My ^ and 
x(0) = 0. It follows that any solution d of the sliding mode 
differential equation monotonically converges to do. At the 
same time, do £ [d_,d+]. Hence d will never leave the 
interval [d_, d + ], the sliding mode will never be terminated, 
the inequality £ + Vt > will never be violated, and d — >■ 
do, d — >■ as £ — >• oo. Application of the Lyapunov's first 
method to the equation y = — x(y) at the equilibrium point 
y = shows that the convergence is exponentially fast. • 

Lemma 6 Both relations d + \{d — do) = 0,£ + Vt > 
become true at some moment following that when the 
distance d reduces to d av . At this moment, d £ [d_, d+] and 
moreover, d > d_ until this moment. 



PROOF. Let t — to be the time when the distance reduces 
to d av . If the required relations hold at t — to, the claim is 
true. Otherwise, d + x(d — do) ¥" f° r t > *o, t ~ to. 

Suppose that d + x{d — do) > for t > t ,t « i - Since the 
control law (6) is ultimately activated by the assumptions 
of Theorem 1, d + x(d — do) arrives at zero at some time 
t\ > to such that d(ti) < d av . By the argument following 
(8), d(t) > do > d- Vt £ [t ,*i], and d(h) £ [d ,eU] C 
[d-, d+\. It remains to note that £ + Vt < is impossible 
by Lemma 4. 

Suppose that d + x(d — do) < for t > to, t « to. Since 
t = to and until the first time t\ when the equation d+x(d— 
do) = becomes true, the vehicle moves with the constant 
control u = —u,u := (V — vq)/L. Now we analyze the 
motion of the robot driven by the constant control u = —u 
for t £ [<o,^o + Ti.5 mm]. By definition, this is a launching 
motion and the vehicle velocity v(t) rotates clockwise at the 
rate u. By the last claim from Assumption 5, there exists 
i* £ [to + T tum,to + r i.5 mm] sucn that as t runs from t 
to t*, the total rotation angle of the vector r(t) — r*(t) 
does not exceed /3 := {t* — T tum — to)u. Meanwhile v(t) 
rotates through the angle 2ir + j3. Hence there exist two time 
instants Si £ [to, to + Ti.s turn] , i = 1,2 such that v(si) and 
( — l)' l [r(si) — r*(si)] are co-linear and identically directed 
fori = 1,2. Thanks to (18), 



(-iyd(si) =v + (-iyV N ( Sl ) >v- \V N \ >' (1 - X v )v 

(13) (16) 

> r] v v > w* > \x[d(si) - do] I 

Thus the continuous function of time S assumes values of 
opposite signs at t = Si,S2. It follows that S inevitably 
arrives at zero within [t , t + Ti.s mm]- Since d £ [d-,d+] 
at this moment by (2) in Assumption 5, Lemma 4 implies 
that (24) cannot be true at this moment. So (23) is true due 
to (25), which and completes the proof. • 



Proof of Theorem 1. This theorem is immediate from Lem- 
mas 5 — 6, along with (22) and (28). The last claim of the 
theorem is straightforward from (22) and (27). 
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